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Abstract 

<N 

A schematic NJL-type model is employed to investigate kaon and pion conden- 
sation in deconfined quark matter in the color-flavor locked (CFL) phase, explicitly 
referring to quark degrees of freedom. To that end we allow for non- vanishing pseu- 
doscalar diquark condensates in addition to the scalar ones which constitute the 
CN ' CFL phase. Color neutrality is ensured by the appropriate choice of color chemical 

^ potentials. The dependence of the free energy in the Goldstone condensed phases 

, on quark masses and charge chemical potentials is found to be in good qualitative 

- in most cases also quantitative - agreement with the predictions obtained within 
the effective Lagrangian approach. 

o . 

1 Introduction 

It is now generally believed that strongly interacting matter at low temperatures and very 
high densities is a color superconductor in the color-flavor locked (CFL) phase [TJ El El H] 
(For reviews on color superconductivity see, e.g., Refs. HH HI IE] ) • For vanishing quark 
masses this phase can be characterized by the equality of three scalar diquark condensates 
in the color and flavor antitriplet channel, 
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S22 = S 55 = S 77 , (1) 

where 

saa' = (q T Cj 5 t a \a< q) ■ (2) 

Here q is a quark field, and ta and A^', A, A' e {2, 5, 7}, denote the antisymmetric Gell- 
Mann matrices acting in flavor space and color space, respectively. In general, these 
condensates are accompanied by induced color-flavor sextet condensates, which are, how- 
ever, small and will be neglected in this article. 

The condensates Eq. (JTJ break the original 577(3) co / or x SU(3)l x SU(3)r symmetry 
of QCD (in the chiral limit) down to a residual SU(3) co i or+ v, corresponding to a common 
("locked") rotation in color and flavor space. As a consequence of the breaking of the 
color symmetry, all eight gluons receive a mass, while the breaking of chiral symmetry 
leads to the emergence of eight pseudoscalar Goldstone bosons. The latter reflect the fact 
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that there is a continuous set of degenerate ground states which can be generated from 
the CFL ansatz, Eq. via axial flavor transformations 

q — > exp(z0 a y7 5 ) q , a = 1, . . . , 8 . (3) 

Under these transformations the scalar diquark condensates, Eq. (J2]), are partially rotated 
into pseudoscalar ones, 

Paa' = (q T C r A X A > q) . (4) 

These condensates are expected to become relevant when we consider the less perfect 
but more realistic situation where chiral symmetry is explicitly broken through non-zero 
quark masses or external charge chemical potentials. In this case the CFL state can 
become unstable against developing non-zero values of Paa', corresponding to pion or 
kaon condensation [§J EH E] . 

This behavior has been predicted within low-energy effective field theories which can 
be constructed systematically for excitations much smaller than the superconducting gap 
A HOI HH H21 HSl HU Basically, this is Chiral Perturbation Theory (%PT) with the 
essential difference to the well-known at low densities that at high densities the 

interaction is weak and the coefficients can be calculated from QCD within High Density 
Effective Field Theory [TO] ITS'] ITS] , which is valid for energies much smaller than the 
chemical potential. 

Let us briefly summarize the leading-order results which we need for later comparison. 
Adopting the notation of Ref . jTT] , the effective meson masses are given by 

M 2 = 2am q m s , = a(m q + m s )m q , (5) 

where we have assumed equal masses for the light quarks, m u = irid =: m q . At asymptotic 
densities, the coefficient is given by a = 3 A 2 j (ir 2 f 2 ) [12] . where A is the CFL gap and 
the pion decay constant is given by f 2 = (21 — 8 ln2)/i 2 / (367r 2 ). 
The mesons experience effective chemical potentials 

2 2 2 2 

m s — m _ m s — m q 
= , &K+ = H ^ — ' = — 2// — ' ^ 

and the same with the opposite signs for 7r~, K~, and K°. The above expressions imply 
that meson condensation takes place if jli exceeds the mass of the corresponding meson. 
In this case the thermodynamic potential of the system is lowered by 



f 2 M 2 

8Q i = -^^(l-cosfl) 2 , cOS# = -rf, (7) 



while for fl 2 < M 2 the CFL state is stable, i.e., 9 = 0. 

Note that the above expressions are of order zero in a s . They are therefore universal 
in the sense that they should not only hold in QCD, but in any model with the same 
symmetry pattern, as long as higher-order corrections in the interaction are small. 

Aim of the present article is to study this mechanism and the corresponding CFL 
+ Goldstone phases within a schematic NJL-type model, explicitly referring to the di- 
quark condensates saa' and paa'- Our main motivation is the fact that NJL-type models 
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have successfully been used to investigate various color superconducting phases and their 
competition (see Ref. [T7] and references therein), but so far the possibility of Goldstone 
boson condensation has not been taken into account*. In particular the recently proposed 
gapless phases at densities below the CFL regime (THl EDI EJj have only been studied 
within NJL-type models which neglect Goldstone condensates. In the present paper we 
want to lay a basis for a more complete analysis in the future. We thereby focus on the 
general mechanism. We investigate whether and under which conditions Goldstone con- 
densation takes place in NJL-type models and compare the results with the predictions 
of the effective theories. To keep the analysis as simple as possible we do not consider qq 
and qij5T a q condensates at the present stage. An extension of the model in this direction 
is straight forward. 



2 Chiral transformations 

Before we define our model we should discuss the axial transformations introduced in 
Eq. (jSJ) in some more details. In this article we consider exclusively flavor non-diagonal 
modes of Eq. (J3J). According to their quantum numbers, these modes can be identified 
with charged pions ("vr" 11 ", a = 1,2), charged kaons ( U K , a = 4,5), or neutral kaons 
("K°", a = 6,7). 

As mentioned earlier, under these transformations the scalar diquark condensates are 
partially rotated into pseudoscalar ones. For a K° mode, for instance, the transformed 
condensates read 

9 

P52 = i sin 2 A ~ ^7) s 22 
9 

p' 25 = i sin - (0 6 + i0 7 ) s 55 

(8) 

where 9 = \/9\ + 9% and 9 a = 9 a /9. The primed and unprimed condensates refer to the 
values in the transformed state and in the original CFL state, respectively. Similarly, 
under K ± transformations, S22 and S77 are partially rotated into P72 and P27, respectively, 
while tt ± transformations, rotate S55 into ^75 and S77 into P57. 

For equal quark masses, because of the residual SU(3) co i or+ v symmetry of the CFL 
state, the same combination of diquark condensates could be reached via axial color 
t r ansf or mat ions 

q — > exp(i6 a -±7 5 ) q , a = l,...,8. (9) 

*For a recent NJL-model analysis of pion and kaon condensation in color non-superconducting quark 
matter, see Ref. [T%] . 



K° : s' 22 = cos - s 2 2 

S77 = S77 
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For instance, for a = 6, 7, this transformation yields 

9 9 
K'° : s' 22 = cos - s 2 2 p'25 = i sin - (0 6 + z6» T ) s 22 

S 55 = COS 2 S 55 P52 = i sin 2 (^6 - «^7) S55 

S77 = S77 (10) 

Hence, as long as Eq. (JTJ) holds exactly, both transformations, Eq. (jSJ) and Eq. (fTUj). lead 
to the same state. On the other hand, in the more realistic case of unequal quark masses, 
there is no exact flavor SU{3)y to begin with and, thus, no exact SU{3) co i or+ y in the 
CFL phase. In this case the scalar diquark condensates S22, S55, and S77 are in general 
not equal and the results of an axial flavor transformation and the corresponding axial 
color transformation will be different. Therefore, in order to distinguish between these 
two transformations, we indicate axial color transformations by a prime, i.e., tt ± , K ± , 
and K'°. 

Looking at Eqs. (jHJ) and (|10j) one might think that a maximally meson condensed state 
corresponds to a transformation angle 9 = ir, i.e., when two of the three scalar diquark 
condensates are rotated away completely, and the two pseudoscalar diquark condensates 
receive their maximum value. This is, however, not true. To see this, it is instructive to 
perform an axial transformation on an idealized vacuum state with (uu) = (dd) = (ss) =: 
<fto- For a K° transformation, as defined above, this state goes over into 

(uu) 1 = 0o, (dd)' = (ss)' = cos^o, (qilsT a q)' = sm9 9 a 4> , a = 6, 7. 

(11) 

Obviously, the condensate with the quantum numbers of a kaon is maximal at 9 = tt/2. 
Therefore, since in a color superconductor quark- ant iquark states and diquark states can 
mix, maximally meson condensed states should correspond to 9 = rr/2 in this well. 
This is consistent with the effective Lagrangian description'''. We will see that it is also 
confirmed by our numerical results. 

3 Model 

We consider an NJL-type Lagrangian with a point-like color-current interaction, 

8 

£ = q(ip-m)q - g ^ (q^Kq) 2 , (12) 

a=l 

where q denotes a quark field with three flavor and three color degrees of freedom, rh = 
diag f(m u , m d , m s ) , and g is a dimensionful coupling constant. 

TThe angle 9 in Eq. (0) corresponds to the rotation angle 7r a // w of the chiral field S = exp(i7r a A a // w ) 
which transforms like S — * LSi?^, where L and R are the SU(3) matrices which transform the left and 
right handed quarks: — > Lq^, qj? —> Rqa- Hence, for L = , the transformation angle for the quarks 
is 9/2, which is consistent with Eq. ©. 
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Performing a Fierz transformation into the particle-particle channel, the interaction 
part can be rewritten as 



2 r 

^11 = -^9^2\{QilbCT A X Al q T )(q T Ci^T A X Al q) + (qCr A X A >q T )(q T Cr A X A/ q) 



+ 



(13) 

where we have listed only those terms which could potentially lead to a condensation in 
the channels s AA i or p AA >. The mean-field thermodynamic potential in the presence of 
these condensates is then given by 

nCT.WJ) = ~TJ2 I 7o^3 \^ ^Q-S-\iu n ,$)) + 2 -g ^T(|W + \paa>\ 2 ) 

n J ^ I A,A> 

(14) 

where T is the temperature, ui n = (2n — 1)ttT are fermionic Matsubara frequencies, and 
{fMfc} (/ £ {«> rf, s}, c G {r, &}) denotes the set of chemical potentials related to the 
conserved flavor and color densities rif c of the Lagrangian. They are often given as linear 
combinations of the chemical potentials /x, /iq, /x 3 , and /i 8 j23], which are related to the 
total quark number density n, the electric charge density uq, and the color densities 
n 3 = n r — n g and = (n r + n g — 2tib)/\/3, respectively (n c := J2f n fc)- 
The inverse propagator in Nambu-Gorkov formalism reads 

o-i,Vi . - f # + A^7°-^ Eaa'( A aa'7s + A%,)^aAa' \ n ,x 

with /t = diagf c (fj,f c ) and the scalar and pseudoscalar diquark gaps 

4 4 
A l4' = A% = — 9 Paa 1 ■ (16) 

In the following we restrict ourselves to the three scalar condensates of the CFL ansatz, 
Eq. ((TJ), plus two pseudoscalar condensates, corresponding to either ti ± (p 57 and 
(p 2 7 and P72), or K° modes (p 25 and p 52 )- In other words, we allow for condensation 
only in one of these modes at a time, which is a reasonable assumption The inverse 
propagator, which is a 72 x 72 matrix, can then be decomposed into a 40 x 40 block 
and two 16 x 16 blocks. Making use of the fact that the integrand in Eq. ()14j) does not 
depend on the direction of the 3-momentum, the dimensionality of the blocks can further 
be reduced by a factor one half. The determinants of the remaining 20 x 20 and 8x8 
matrices and the Matsubara sum are evaluated numerically. 

Until this point, the thermodynamic potential depends on our choice of the diquark 
gaps. As standard, the stable selfconsistent solutions are given by the values of A S AA , and 
A P AA , which minimize Q. In most cases below, however, we will restrict the minimization 
procedure to the subspace of condensates which can be obtained from the CFL solution 
(A^, = 0) via chiral rotations, Eq. (JSJ) or Eq. (J0J), in order to compare our results with 
the effective Lagrangian approach. 
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4 Numerical results 



In most of our numerical studies we adopt the model parameters of Ref. |22|, namely 
a 3-momentum cut-off A = 600 MeV and gA 2 = 2.6 for the coupling constant. For 
comparison with the predictions of the effective Lagrangian approach we keep the quark 
masses as variable parameters. As before, we assume equal masses for up and down 
quarks, m u = m& =: m q . All calculations are performed at zero temperature and at fixed 
H = 400 MeV. 

To have a well-defined starting point, we begin our analysis in the chiral limit, m q = 
m s = and with equal chemical potentials for all quarks, i.e., \iq = fi 3 = fi 8 = 0. For this 
case we find a CFL solution with A 22 = Ag 5 = A 77 = 104.2 MeV. As expected, this is not 
a unique solution for the ground state, but an infinite set of degenerate ground states can 
be obtained via chiral rotations. In particular, any rotations in n^, K^, or K° direction 
leave the ground state free energy invariant. 

Next we introduce a finite strange quark mass m s = 120 MeV, while leaving the up 
and down quarks massless. The explicit breaking of the SU(3) symmetry causes a slight 
asymmetry in the ground state, i.e., instead of an ideal CFL solution, Eq. ((TJ), we get a 
5% splitting of the gap parameters: A s 22 = 105.4 MeV and A s 55 = A s 77 = 100.7 MeV. It 
turns out that the corresponding free energy is still invariant under chiral rotations into 
the ir^ direction, whereas it becomes disfavored, if we perform a rotation into the kaonic 
directions. At first sight, this result seems to contradict the effective Lagrangian results, 
as summarized in Eqs. © - (J7J). 

At this point we should notice that, by construction, the effective Lagrangian contains 
only colorless states (as a result of integrating out static gluons, see, e.g., Ref. [13]) whereas 
the above ground state of the NJL model is not color neutral. In fact, it is well known 
that the introduction of unequal masses in the CFL phase leads to colored solutions in 
NJL-type models, unless this is corrected for by introducing appropriate color chemical 
potentials [23] • This can be thought of as an effective way of simulating static gluon 
background fields. In the above example we need to introduce a chemical potential /ig = 
—10.12 MeV in order to obtain a color neutral CFL solution. This solution has a smaller 
splitting of the gaps, A s 22 = 103.4 MeV and A s 55 = A s 77 = 101.8 MeV, and is about 
0.2 MeV/fm 3 higher in free energy as the colored one. 

However, the essential observation is that, restricting ourselves to color neutral solu- 
tions, the CFL state does no longer correspond to the ground state of the system if we 
allow for chiral rotations. This is illustrated in the left panel of Fig. ^ where the relative 
change of the free energy density is plotted as a function of the angle 9 for chiral rotations. 
While the rotations in pion direction (dashed line) still leave the free energy invariant, the 
kaonic transformations now lead to a reduction of Q. (Note that charged and neutral kaon 
solutions are degenerate under the present conditions.) Here the dotted line indicates the 
result of an axial flavor rotation (K ± or K°), and the solid line the result of an axial color 
rotation (K or K °). As one can see, the latter is more favored and symmetric about a 
minimum at 9 = 7r/2, while the former is less favored and slightly asymmetric. 

This can be explained by the fact that, e.g., a K° rotation transforms the larger S22, 
in which u and d quarks are paired, into the us condensate P52, whereas the smaller 
us condensate s 55 is rotated into the ud condensate p 25 . Hence, the number of strange 
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Figure 1: Relative change of the free energy density as a function of the chiral angle for 
m u = = and m s = 120 MeV. The calculations have been performed at u = 400 MeV 
and [iq = (left) or fiQ = — 50 MeV (right). The color chemical potentials have been adjusted 
to obtain color neutral solutions. The various lines correspond to different modes as indicated 
in the figure. The crosses indicate the free energy obtained by minimizing the thermodynamic 
potential for color neutral matter without restriction to chiral transformations. 



quarks is larger at = n than at = 0, which explains why = tx is less favored. In 
a K transformation, on the other hand, ud condensates remain ud condensates and ds 
condensates remain ds condensates, because axial color transformations do not change 
the flavor structure. Therefore the number of strange quarks at = ir is equal to the 
number of strange quarks at = 0, and both states are degenerate. 

None of the above results corresponds to the real minimum of the thermodynamic 
potential, because we have restricted ourselves to certain chiral rotations of the CFL 
solution without allowing for "radial" variations. This means, instead of freely varying 
the scalar and pseudoscalar condensates under consideration, we have linked them to the 
CFL solution and a single parameter 0, e.g., as given in Eqs. (JBJ) and (fTUJl for the K° or 
the K case, respectively. If we abandon this constraint we find solutions with even lower 
free energies, indicated by the cross in the figure. These solutions are very similar to a K 1 
state at = ir/2, in the sense that, e.g., |A| 2 | = (A^l and |A§ 5 | = (A^l, but the values 
of the ud gaps are slightly larger (74.2 MeV = 104.9 MeV/v^2) and the values of the us 
gaps are slightly smaller (70.9 MeV = 100.2 MeV/\/2) than the rotated CFL gaps. 

In fact, in our context the .^'-transformations should also be interpreted as certain 
"radial" transformations on top of the axial flavor transformation K, and not as axial 
color transformations. Note that, even in the chiral limit, axial color transformations, 
Eq. 0, are neither a symmetry of QCD nor of our NJL- model Lagrangian, Eq. ()12j) . 
They are just a symmetry of the mean-field thermodynamic potential, as long as we 
restrict ourselves to the given set of diquark condensates. It will no longer be the case if, 
e.g., qq condensates are taken into account. In the present model, the difference between 
the K and K' transformations is due to the non-equality of the diquark condensates in the 
CFL phase (cf. Eqs. (JHJ) and (jlOj)). This is a higher-order effect, which is not taken into 
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account in the leading-order results of the effective Lagrangian approach. Hence K and 
K' transformations are in principle equally good starting points for a comparison with the 
effective theory. However, in some aspects the K' transformations behave more similar to 
the kaon modes in the effective theory since they are symmetric about 6 = 7r/2*. Moreover, 
they come closer to the true minimum of the thermodynamic potential. Therefore, for 
most of quantitative studies in this article we perform K' transformations rather than K 
transformations. Note that for the pion modes there is no difference anyway. 

In the right panel of Fig. ^ we show the result of a similar analysis, again for m s = 
120 MeV, but now for a non-vanishing electric charge chemical potential fiQ = —50 MeV. 
The color chemical potentials 113 and 11% have again been adjusted to ensure color neutrality 
at each point. Under these conditions all flavored Goldstone modes lead to a reduction 
of the free energy, but the pionic mode (dashed line) is the most favored one, followed by 
the charged kaons (dash-dotted). The neutral kaon is not sensitive to jiQ and shows the 
same behavior as in the left panel. As before, none of the modes corresponds to the real 
minimum of the thermodynamic potential one gets if the condensates are not constrained 
to certain chiral rotations. This free energy is again indicated by a cross in the figure. 

At least qualitatively, the results presented in Fig. ^ are in good agreement with the 
predictions of the effective Lagrangian approach. For a more quantitative comparison 
we again restrict ourselves to chiral transformations of the CFL solution. In Fig. |2] we 
show the behavior of the free energy gain of the meson condensed solutions (6 = 7r/2) 
for massless up and down quarks as functions of \xq (left panel) and m s (right panel). 
At each point, the color chemical potentials /X3 and /Xg have been tuned to ensure color 
neutrality. 

The left panel shows the behavior of SQ for the three Goldstone boson condensates 
ir^, K , and K at fixed chiral angle 9 = 7r/2 as functions of \xq at fixed m s = 120 MeV. 
Since, according to Eqs. (J5J) to (J7J), we expect 5Q n to behave like /j,q, we plot the di- 
mensionless ratio <5f2/(/i 2 /ig). Indeed, our result for the pionic mode is constant to a 
very high degree (dashed line). 5Q K ' , on the other hand, does not - and should not - 
depend on /j/q, and hence the corresponding curve (solid line) decreases like 1/7-tg in the 
figure. Finally, the charged kaons behave more complicated. According to Eq. (JSJ), \i K ± 
vanishes at \iq = — m s /(2/x), leading to a maximum in the free energy. Qualitatively, this 
is well reproduced in our model calculation (dash-dotted) line. In fact, in the \xq interval 
between -13 MeV and -22 MeV, 5Q K >± even gets slightly positive. This seems to indicate 
that the kaon is not exactly massless but has a small mass of about 5 MeV. Then, since 
the leading-order (a°) effective theory predicts m K = for m q = (see Eq. (JHJ) ) , this 
must be a higher-order effect due to interactions. There are, for instance, QCD correc- 
tions proportional to a s m 2 s [213, which have, however, the wrong sign. Nevertheless, there 
could be other terms which lead to the observed behavior. We will come back to this 
issue below. 

For comparison, we also show the effective Lagrangian results for the free energy as 
given by Eqs. (JSJ) to (JZJ) (dotted lines). Here we treat f n as a free parameter which 

•fFor instance, the color chemical potentials needed to neutralize the K' solution at the minimum at 
8 = t: / 2 satisfy the relation = v^Ms- This is in nice agreement with the ratio between the third and 
the eighth color component of the static gluon field, = \fi 4% , obtained in Ref. 24_ for neutral CFL + 
K° solutions. This relation does not hold for the K solution in our model. 
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Figure 2: Free energy density difference between meson condensed phase (8 = tt/2) and CFL 
phase (Q = 0) as functions \xq (left) or m s (right). The calculations have been performed for 
m u = rrid = and at fixed [i = 400 MeV. The color chemical potentials have been adjusted to 
obtain color neutral solutions. In the left panel m s = 120 MeV, in the right panel [1q = 0. The 
various lines correspond to different modes as indicated in the figure. The dashed curve labeled 
"K'(1.5)" in the right figure has been calculated with a reduced coupling constant gls? = 1.5. The 
bold lines have been obtained within the present model while the thin dotted lines correspond 
to Eqs. © to © with f n = 80 MeV. 



is fitted to reproduce the pion. We obtain f n = 80 MeV, which also gives reasonable 
fits to the kaons. The fitted value agrees quite well with the leading-order result f n = 
y/(21 - 81n2)/V(67r) = 83.4 MeV. A gain, there should be corrections to this formula due 
to the interaction. These could depend on the meson channel as well as on m s and \Lq. Of 
course, meson masses and decay constants could be calculated explicitly within the NJL 
model, summing up qq loops in Nambu-Gorkov formalism (which are essentially diquark 
loops) and coupling them to an external axial current. This is left for future work. 

In the right panel of Fig. |2]the free energy of the kaonic modes (K and K') at fiQ = 
is shown as a function of m s . For comparison we show again the result corresponding to 
Eqs. (0) to (0) for f n = 80 MeV (dotted line). Since according to these equations we expect 
5Qk to behave like m\ we plot the ratio SQ/m^. However, whereas for m s > 100 MeV 
this ratio indeed is roughly constant, we find strong deviations at lower strange quark 
masses. In particular, below m s = 31 MeV the kaon modes are no longer favored. 

Again, this indicates the presence of higher-order interaction effects. Note that the 
deviations from the m\ behavior show up in a regime where m s is smaller than the CFL 
gap A. Hence, terms of the order m^A 2 , A 4 or m 2 A 4 //x 2 j2H], which are exponentially 
suppressed in the weak-coupling regime, could be parametrically comparable to m\. To 
check this point, we redo the calculation with a reduced coupling constant, gA 2 = 1.5. 
The CFL gaps are then about 30 MeV, instead of ~ 100 MeV which we had before. The 
resulting Sfl for the K' mode is indicated by the dashed line in the figure. Obviously, 
this curve shows an almost perfect m 4 behavior down to m s = 25 MeV and there is 
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Figure 3: Chiral angle (upper panels) and corresponding free energy gain (lower panels) of 
meson condensed phases at \x = 400 MeV for non- vanishing light quark masses m q . The full 
squares, circles, and triangles indicate the NJL model results for m q = 5 MeV, 10 MeV, and 
15 MeV, respectively. The dotted lines are the corresponding results of Eqs. (JSJ) to (J7J) with 
a = 0.169 and f n = 80 MeV. Left: pion condensate for m s = 120 MeV and varying fiQ. Right: 
kaon condensate (K 1 ) for [xq = and varying m s . 



no indication that kaon condensates become disfavored at small m s . (Unfortunately, we 
cannot go to even smaller values of m s because the numerical uncertainties become too 
large.) This means, there is at least numerical evidence that the effect seen with our 
original parameters is indeed due to interactions. 

Finally, we consider non-vanishing masses for both, strange and non-strange quarks. 
Then, according to Eq. (jSj), the effective Goldstone masses become finite even at leading 
order, and m s or fig have to exceed certain threshold values to enable Goldstone boson 
condensation. Our results are summarized in Fig. El The full squares, circles and tri- 
angles correspond to NJL-model calculations with m q = 5 MeV, 10 MeV, and 15 MeV, 
respectively. 
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In the two left figures we show the behavior of the pionic solutions for fixed m s = 
120 MeV as functions of \xq. In the upper panel the chiral angle 9 is shown which minimizes 
the free energy. As expected, this angle is zero for small values of \hq\ and becomes non- 
zero above a threshold value which depends on the quark mass. This behavior can be 
described well by Eqs. © to (0) if we treat the constant a which enters into the expression 
for the meson masses as a free parameter. The dotted lines, which run almost perfectly 
through the points, have been obtained with a = 0.169. However, this is only about 1/3 of 
the leading-order result a = 3A 2 /(n 2 f 2 ) Q2|, if we take A ~ 100 MeV and f v = 80 MeV. 
Again, this discrepancy must be an interaction effect. Nevertheless, since in Fig. EJwe 
found the leading-order results to be in rather good agreement with the pion channel, 
a factor 3 in the present case is a bit surprising. To get some insight, we have again 
repeated the calculation with gA 2 = 1.5, i.e., A ~ 30 MeV. For this case we find the 
fitted value for a to be about 2/3 of the leading-order prediction. This seems to indicate 
that the latter is reached asymptotically, but the convergence is relatively slow. A more 
systematic investigation of this point is certainly necessary. 

In the lower left panel we display the free energy gain (for gA 2 = 2.6). The dotted lines 
indicate the results of Eqs. (0) to (JJJ), taking the previously fitted constants a = 0.169 
and f n = 80 MeV. Obviously, we obtain a perfect description of the NJL model results 
without refitting these constants. 

In the right two figures the analogous quantities are shown for kaon (K') condensates 
at jiQ = as functions of the strange quark mass. The dotted lines correspond to Eqs. (JHJ) 
to (0) with the values for a and f n fitted in the pion channel. Obviously, the overall 
behavior is well reproduced, although the quantitative agreement is not as good as in the 
left two panels. In fact, the deviations are consistent with our earlier results: As one can 
see in the upper figure, the NJL-model results for the thresholds for kaon condensation 
are at somewhat higher values of m s than the dotted lines suggest. This could again 
be explained by a positive correction term to the kaon mass beyond the leading order. 
Accordingly, the NJL-model results for Sfl (lower right panel) are also shifted to higher 
values of m s as compared with the dotted lines. On the other hand, for large m s , the light 
quark masses become less relevant and the results approach the m q = values of Fig. |21 

5 Summary and outlook 

We have studied pion and kaon condensation in the CFL phase within an NJL-type 
model. To that end we have performed a mean-field calculation allowing for non-vanishing 
expectation values of certain scalar and pseudoscalar diquark condensates. Main focus of 
the present article was on the general principle. We have explicitly shown that Goldstone 
condensates, i.e., non- vanishing expectation values of pseudoscalar diquark condensates, 
develop in reaction to a finite strange quark mass or a finite electric charge chemical 
potential. In this context it was essential to introduce color chemical potentials to ensure 
color neutrality of the CFL and the CFL + Goldstone solutions. 

More quantitatively, we found good over-all agreement of the NJL-model results with 
predictions obtained within chiral perturbation theory in the CFL phase to leading order 
in the interaction. The most intriguing exception is the coefficient a, which determines the 
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dependence of the Goldstone boson mass on the quark masses. Fitting this coefficient to 
our results, it turned out to be only one third of the leading order High Density Effective 
Theory value if the coupling is relatively strong (A ~ 100 MeV) and two thirds for 
a weaker coupling (A ~ 30 MeV). The leading-order value may thus be reached at very 
weak couplings only, but this deserves further investigations. Apart from this point we 
found some deviations in the kaon sector which we also identified as interaction corrections 
to the leading order. 

Having demonstrated the general consistency of the method, we should now extend our 
analysis to perform a more complete study of the phase structure of strongly interacting 
matter at high densities. In the present article, we have in most cases restricted the 
space of possible diquark condensates to certain chiral rotations, in order to facilitate the 
comparison with the effective theories. This restriction should be relaxed. As shown in 
Fig. ^ this can lead to a further reduction of the free energy. 

The model should also be extended to include quark-antiquark condensates, like (uu), 
(dd), and (ss), which lead to density dependent effective quark masses and in this way 
could influence the phase structure considerably [221 • Now, in order to study Goldstone 
boson condensation, we also need to take into account the corresponding chiral partners, 
(qi^Taq). These condensates are the essential degrees of freedom to describe pion and 
kaon condensates in color non-superconducting phases, which have recently been studied 
in great detail in Ref. [IB]. In the CFL phase, the dominant effects should be due to the 
diquark condensates, but nevertheless, quark-antiquark condensates could lead to impor- 
tant corrections. We should also note that our toy- model Lagrangian, Eq. (|12|). misses 
instanton effects, which could significantly contribute to the Goldstone boson masses j26j . 
Finally, it would be interesting to extend the analysis to the gapless CFL phase. 
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